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Abstract. - The reversion of the time evolution of a quantum state can be achieved by 
changing the sign of the Hamiltonian as in the polarization echo experiment in NMR. In this 
work we describe an alternative mechanism inspired by the acoustic time reversal mirror. By 
solving the inverse time problem in a discrete space we develop a new procedure, the perfect 
inverse filter. It achieves the exact time reversion in a given region by reinjecting a prescribed 
wave function at its periphery. 



Performing a time reversal experiment in a quantum system seemed an imposible task until 
Erwin Hahn realized that this would "only" require inverting the sign of the Hamiltonian. This 
was achieved for a set of independent nuclear spins precessing in a magnetic field through the 
sudden application of a radio frequency (rf) pulse and became the basis of the spin echo [3]. 
When a similar strategy was applied to a system of interacting spins, a true many-body 
system [4], the first practical realization of a Loschmidt daemon [1,2] was finally achieved [5]. 
We dub this procedure a hasty daemon, as it involves the global and instantaneous action 
of a rf pulse. Notably, this action tends to be quite dependent on the underlying instability 
(chaos) of the corresponding classical system [6]. 

On the other hand, during the last decade, Fink and his group developed an experimental 
technique, called time-reversal mirror (TRM), which reverse the propagation of acoustic waves 
[7]. A pulse at a point source ro on a working region C is detected as it arrives to an array 
of transducers at positions r^, typically surrounding the cavity C. Their registries ip(ri,t) 
are recorded until a time tji at which the amplitudes ip(ri,t) have become negligible. These 
transducers can act alternatively as microphones or loudspeakers. Afterwards, each one re- 
emits in the time reversed sequence i.e. producing an extra signal Xinj. ( r ii tn+5t) = cipfa, tji — 
St) where c is controlled by the "volume" knob. The experiments show that these waves 
tend to refocus at the source point at the time t = 2tu, i.e. a Loschmidt echo is formed! 
The robustness of the time-reversal procedure is a prominent feature of these experiments. 
Surprisingly, systems with a random scattering mechanism are particularly stable. In fact, 
there is a precise prescription for time reversal [8] that requires the control of the field function 
and its normal derivative over a surrounding surface. However, the TRM procedure is quite 
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effective even when these requirement are not completely satisfyied. On this basis, several 
applications in communications [9] and medicine [10] were already performed. It is clear that 
these experiments introduce a different procedure for time reversal: a persistent action at the 
periphery that we call a stubborn daemon. At this point the the first question to ask is: Can 
this concept be applied to a Quantum Mechanical system? 

A quantum experiment would touch an issue usually overlooked: the use of the Schrddinger 
equation (SE) with a time dependent source. This phenomenon is not merely academic since 
it appears in many areas: the gradual injection of coherent polarization in the system of 
abundant nuclei through an NMR cross-polarization transfer [11]; the creation of a coherent 
excited state [12] through particular sequences of laser pulses at slow rates of pumping; and 
an a.c. electrical conductivity experiment where the electrodes are fluctuating sources of 
waves [13]. However, there is no general answer [14,15] to the "inverse time problem": What 
wave function must be injected to obtain a desired output? In what follows we solve this 
problem for a reasonably general case and use that solution to implement a protocol for a 
perfect quantum time reversal experiment. 

Let us consider a one dimensional system with a wave packet tp(x,t) localized around 
a point x at time t = and traveling rightward. The first question to answer is: If we 
record the wave function as a function of time only at a point x s > x , is it possible to use 
this information to recover the same wave function? Answers to this question were given in 
Ref. [15] for some particular potentials. To solve this problem in a general way, we resort to 
a discrete Hamiltonian: 

« = E EjSfa + E^.,-,^-^ + 

3 3 

here c^ and £j are the creation and annihilation operators for a particle at the coordinate 
Xj = ja where a is the lattice constant. In the standard notation [?] the kinetic energy yields 
the hopping term V x . x . = —V8i±ij with V — h 2 / {2ma 2 ), while the potential energy U(x) 
fixes the "site energy" Ej = U(xj) + 2V. We want to express ip(x = ja,t) = (0\dj \ip(t)) in 
terms of the wave function at the position x s of the detect or /source. We start with the usual 
expression: 

ii{x,t) =J2^G* Xn {t - t Mx n ,t ), (2) 

n 

where the time retarded Green's function G XXn (t— to) satisfies ih-^G XXn (t)—^2. H x , Xi G x . Xn (t) — 
dx,x n 5 [t\- In the energy representation: 

i>(x n ,t ). (3) 

At this point, we separate the space in two portions: one will be the working space C where 
one intends to control the wave function. The other, the outer region, is the complementary 
infinite region that contains the scattering states. Note that in our discrete calculation the 
connection between both regions is achieved by the hopping term V x Xs , connecting the 
sites x s and x s+ i at both sides of the boundary. We use the Dyson equation G X Xn = G X Xn + 
Gx,x s+1 Vx s+1 ,x s G XstXn relating the Green's functions G x . x ., of the semi-spaces defined by 
V Xs+i Xb = 0, with the complete one. Hence, for x > x s 

^{x,t)=ihY, I ^exp[^^]{G« Xi+1 ( e )K. +1>x J[G*, a;n ( e )V'(a: n ,to)]. (4) 

n J 
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The sum within square brackets can be identified with the energy representation of the wave 
function ( i.e. t/j(x s ,s)= J2 n [Gx s ,x n ( £ )V'( a; Ti, to)]). Besides, when x > x s the Dyson equation 

becomes G XXs (s)= (g xXs+1 (e)V Xs+iX \G XsXb . From this, we evaluate the term in curly 
brackets that we replace in the Eq. J3J to obtain, for x > x s : 

ip(x,t) = ifi J 5^-exp 

The SE with a source, ih-^ip(x, t) — J2 X 'M XiXn ip(x n ,t)= Xinj.(x, t), has the general solution: 
i,{x,t)=\K f ^exp[^]^G« X3 ( £ ) Xi „j.(x s , £ ). (6) 

x B 

This allows us to identify 

X S ource( X s> £ ) = QR 1 ^ 1p(x s e), (7) 

as the Fourier transform (FT) of the function that must be injected at each instant in order to 
obtain the target function. This result is valid in any dimension and for an arbitrary potential 
U(r). The condition x > x s becomes r G C, and one must interpret ip(x s ,e) as a vector 
whose components are the wave amplitudes at the N sites {r s } defining the boundary B . 
Similarly, one recognizes G x x (e) as the N x N matrix providing the correlations between 
these sites. To our knowledge this is the first solution to the inverse time problem. The key 
feature allowing this simple solution was the representation of the Schrodinger equation in 
a discrete basis. This enabled a natural separation into complementary subspaces that are 
re-connected through the Dyson equation. 

Time-reversal via injection. In the following, we propose a gedanken scheme to achieve 
a perfect time-reversal of an arbitrary wave packet by assuming that a persistent non invasive 
injection and detection of waves at a single point is possible. In such conditions one would 
create an efficient stubborn daemon: the Perfect Inverse Filter (PIF). We illustrate this by 
considering an incoming wave packet in a semi-infinite space bounded by an infinite barrier at 
imax. which, together with a scattering barrier, define a reverberant region ( see Fig. [TJ. At 
the point x Sl located to the left of the scattering barrier, we alternate the use of an injector 
and a detector of wave function (probability and phase). This set up is a particular realization 
of the "sound Bazooka" scheme implemented by Fink's group. However, instead of using the 
TRM, we proceed as follows: 

1) We calculate the response of the system to an instantaneous excitation at site x s i.e. 
G XstX3 (t) and compute its FT, G^ Xs {e). 

2) We start with the empty cavity ( %j)(x > x s ,0) = 0) and a wave packet that travels 
towards it (e.g. a Gaussian centered at xq < x s ). The probability density at time zero is 
shown in the top of the left panel of Fig. [3 It is followed by a sequence of snapshots of the 
density at selected times in the range [0, 2£r] increasing from top to bottom and continuing in 
the right panel from bottom to top. The inject ion/ detect ion point is indicated by a vertical 
dotted line in each panel. 

3) During the period (0,T™ = t R - ti) the wave packet performs a free evolution: it 
enters to the cavity, collides with the barrier and then bounces back in the wall at the right 
end of the system and finally escapes towards the outer region at the left side. See left panel 
of Fig. Provided that there are no localized states in the cavity and that the wave packet 
that escapes to the outer region does not return, the condition ip(x s ,t > tu) ~ can be 



Gl Xa {e) ( M x s e). (5) 
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fulfilled. During the whole period T^Jf = t R — 1\ the wave function amplitude and phase at 
x s are registered. The range (ti, tu) should contain the support of iJj(x s , t), i.e. ip(x. t) = for 
x > x s and t £ (ii,f#). 

4) Now our target function is tp rev (x s ,t R + St) — ip*(x s ,t R — St) with < St < T*^f, 
i.e. the wave packet with reversed evolution. Using the information registered in the previous 
step, we Fourier transform it 



and normalize it according to Eq. JJJ. Transforming back to time we get the actual time 
dependent injection acting for a time T^f '. The injection also produces a wave packet that 
travels to the left, i.e. escaping to the left outer region, see Fig. |3] Hence, perfect time 
reversion is restricted to the cavity, i.e. x > x s . 

5) After injection has seased, the original wave-packet is recovered at time 2t R with an 
inverted momentum: this is the Loschmidt Echo. Figure [21 also shows, in dotted line, the 
echo resulting from TRM procedure [7], which in this case would require the recording only 
the outgoing wave described in step 3) which is time-reverted and reinjected without further 
procesing. 

In Fig. the density at site x s is shown at different times. The actual PIF density 
is plotted with a solid line while the injected density at each time is shown with a dashed 
line. Notice that the injection intensity has to provide a wave propagating toward both the 
cavity and outer region. We also show, with a dotted line, the density obtained from the 
TRM procedure. Note that this density also exhibits an Echo at time 2t R but with a reduced 
amplitude as compared with the original signal. While the PIF and TRM densities differ in 
their magnitude, their shape in [tn, 2tp\ is remarkably simmilar. This indicates the stability of 
the TRM in the condition considered. We emphasize that by injecting probability amplitude 
according to Eq.Q we exactly reverse the forward evolution of the initial wave packet. The 
correction imposed by the PIF procedure becomes non-trivial in cases where the incoming and 
outgoing signals superpose. In these cases the PIF procedure "filters" the outgoing portion 
as can be appreciated in Fig. In the upper panel we display the forward evolution as 
registered at site x s . In this case, one can roughly indentify three time regimes: entrance 
(IN), escape (OUT) and a Mixed region when both components interfere. The lower panel 
displays the time reversal procedure. The TRM would inject only the OUT portion of the 
registry shown in the upper panel. In contrast, the PIF procedure yields an injection that 
extends to the mixed region shown ligth grey shaded (yellow area on-line). The dark shaded 
(blue area on-line) PIF intensity constitutes a substantial improvement over the gray shaded 
(cyan area on-line) TRM signal. 

The PIF protocol is valid for the reversal of any scalar waves as long as they satisfy a 
linear equation. Then, different propagators are described by the Green's functions. The basic 
ingredients apply to elastic or electromagnetic waves [16] extending the range of applicability 
of the concepts introduced here. 

The implementation of a stubborn Loschmidt daemon in a quantum system is not a simple 
task. However, standard pulsed NMR has the tools. In an ensemble of linear molecules, the 
interactions between nuclear spins can be manipulated [18] to obtain a polarization which is the 
square modulus of a single particle wave function or polarization amplitude [19] and constitutes 
a pseudo-pure state [17]. Detection at each time involves an ensemble measurement and a new 
experiment. In order to generate a local source/detector one resorts to the interaction between 
different nuclei, which can be engineered at will, e.g. a 13 C acts as such probe for a labeled 1 H. 
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In fact, we have been able to inject a wave packet in a 1 H ring and follow its dynamics detecting 
simultanously the amplitude and relative phase at the labeled 1 H [20]. This is a double-slit 
like experiment that allows interferometry in the time domain [21]. If most of the polarization 
stays in the 13 C, it is the small portion transfered to the proton system the one described by 
the theory above. While a full implementation of the quantum TRM or PIF requires setting 
many important experimental details, every step towards that goal would have potential use 
in spectral edition and quantum information processing. More immediately, classical wave 
systems, could benefit from our procedures which can be incorporated in a straightforward 
manner. 

In summary, we have studied the Schrodinger equation with source boundary conditions. 
We have obtained a general solution for the inverse time problem which is expressed in terms 
of the Green's function at the boundary region. Our results enabled us to develop the Perfect 
Inverse Filter protocol to implement a stubborn Loschmidt daemon. It allows one to achieve 
the perfect time reversal of the wave dynamics obtaining a Loschmidt Echo. In some cases, 
this protocol could improve the experimental procedure implemented with sound waves. Now 
that a perfect reversion can be obtained, a number of questions relevant to the Quantum 
Chaos field become pertinent related to the assestment of infidelity sources. On the view of 
the experimental results in sound waves [7], a stubborn daemon yields more robust results 
than its hasty counterpart. Hence, a whole field of study opens up. 
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Fig. 2 - Distribution of probability density (in the cavity and outer region) for different times. Panel 
a) shows the forward evolution between and tR. The backward evolution is shown in panel b). 
The solid line is the result of an appropriate injection, Eq. (J, while the dotted line is obtained 
by injecting only the time-reversed wave recorded at x a during the forward evolution. The initial 
gaussian wave packet, which is centered at x s — 200a, has a/a — 50 and koa — 1. 
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Fig. 3 - Probability density at site x a as a function of time. The actual density (solid line) results 
from the injection of the density shown with a dashed line corresponding with the FT of Eq.Q. The 
density corresponding to the injection of the time-reversal of the recorded amplitude is also shown 
with a dotted line. The injection and recording periods, at site a; s , are specified in the top of the 
figure. 
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Fig. 4 - Upper panel: density at x s in a forward evolution of a wave packet. Lower panel: At rigth, 
the density injected at x s by the PIF procedure (shaded light-grey /yellow on-line) . Its evolution 
determines the correct reversed density (shaded dark/blue on line). The density obtained from TRM 
is shaded grey/cyan on-line. System parameters: Barrier height 0.2V between sites x s + 100a and 
x 3 + 105a and a; max . = 200a. Gaussian wave packet centered at x s — 100a, with a /a — 3 and koa = 0.8. 



